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ABSTRACT 

We present HI synthesis maps of the edge-on starburst NGC 5433 and its environ- 
ment, obtained with the VLA in its C and D configurations. The observations and 
spectral model residuals of the main disc emission in NGC 5433 reveal 3 extraplanar 
features. We associate 2 of these features with coherent extraplanar extensions across 
multiple spectral channels in our data, including a complete loop in position-velocity 
space. Interpreting the latter as an expanding shell we derive a corresponding input 
energy of 2 x 10 54 ergs, comparable to that for the largest supershclls found in the 
Galaxy and those in other edge-on systems. NGC 5433 is in a richer environment 
than previously thought. We confirm that KUG 1359+326 is a physical companion 
to NGC 5433 and find two new faint companions, both with Minnesota Automated 
Plate Scanner identifications, that we label SIS-1 and SIS-2. Including the more dis- 
tant IC 4357, NGC 5433 is the dominant member of a group of at least 5 galaxies, 
spanning over 750 kpc in a filamentary structure. A variety of evidence suggests that 
interactions are occurring in this group. While a number of underlying mechanisms 
are consistent with the morphology of the high-latitude features in NGC 5433, we 
argue that environmental effects may play a role in their generation. 
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1 INTRODUCTION 

High-latitude extensions from the interstellar medium (ISM) 
in spiral galaxies are important probes of disc stability and 
energetics, halo dynamics, and the physical processes that 
drive galaxy evolution. These features are best detected 
in edge-on galaxies, in which extraplanar emission can be 
clearly separated from the underlying disc. Discrete fea- 
tures (e.g. arcs, plumes, and filaments in various ISM trac- 
ers, or 'supershells' in H I as well as broader-scale structures 
such as thick discs or haloes) have now been identified in 
many galaxies using this approach. This has led to ideas 
of galactic circulation involving 'fountains', 'chimneys', or 
possibly outflowing winds which may enrich the intergalac- 
tic medium (Shapiro and Field 1976, Bregman 1980, Nor- 
man and Ikeuchi 1989, Heckman, Armus & Miley 1990, Bre- 
itschwerdt, Volk & McKenzie 1991, Breitschwerdt, McKen- 
zie & Volk 1993). 

Disc-halo features are usually attributed to underlying 
drivers such as supernovae and stellar winds. However, these 
processes alone may be insufficient to produce the highly en- 
ergetic HI supershells seen in some galaxies, or to account 



for the absence of detectable star formation regions or rem- 
nants in others (e.g. Rhode et al. 1999, Perna & Gaensler 
2004). Thus, while star formation may be a key ingredient 
in the disc-halo phenomenon, it may not, alone, be suffi- 
cient to drive disc-halo dynamics. Impacting clouds have 
long been invoked as a possible alternative (Tenorio-Tagle 
& Bodenheimer 1988, Santillan et al. 1999), but galaxies 
that are apparently isolated also sometimes display large 
supershells (King & Irwin 1997, Lee & Irwin 1997). This 
has led to speculation about other possible contributors, in- 
cluding gamma-ray bursts (Efremov, Elmegreen & Hodge 
1998, Loeb & Perna 1998), inflated bubbles from radio jets 
(Gopal- Krishna & Irwin 2000) and magnetic fields (Kamaya 
et al. 1996). Environmental effects have also been invoked 
to explain a variety of halo phenomena, such as differing ex- 
tents of radio (Dahlem, Lisenfeld & Golla 1995) and X-ray 
haloes (Wang, Chaves & Irwin 2003) among various edge-on 
systems. 

This paper represents an attempt to further character- 
ize the importance of environment in driving disc-halo dy- 
namics. A first problem is the lack of information about 



© 2004 RAS 



2 Spekkens, Irwin and Saikia 



the environments of galaxies showing disc-halo outflows. 
We have therefore obtained maps of the large-scale HI in 
the edge-on galaxy NGC 5433 using the Very Large Ar- 
ray (VLA) in its C and D configurations. We examine the 
main disc and extraplanar H I in NGC 5433 using the higher 
resolution C configuration data, and combine the C and D 
datasets to probe its environment. 

NGC 5433 (D = 65 Mpc; H = 70 kms^Mpc' 1 ) is 
an infrared bright spiral (Soifer et al. 1989) that was ob- 
served in the radio continuum surveys of Irwin, English & 
Sorathia (1999; see their fig. 19) and Irwin, Saikia & English 
(2000; see their fig. 1). The 20 cm images from these stud- 
ies reveal a thick disc which resolves into discrete disc-halo 
features at high resolution. Single dish HI observations of 
NGC 5433 have been obtained by Mirabel & Sanders (1988) 
and Schneider et al. (1990) and a low sensitivity VLA inte- 
grated H I map has been presented by Thomas et al. (2002) . 

An optical DSS image of NGC 5433 and its southwest 
(SW) environment is shown in Fig. NGC 5433 appears 
slightly curved, with an outer ring which is distorted on the 
north side. Its RC3 (de Vaucouleurs et al. 1991) classifica- 
tion as an Sdm has been revised to an earlier type by Dale et 
al. (2000), who include NGC 5433 in an Infrared Space Ob- 
servatory Key Project study of normal star-forming galaxies 
(see also Malhotra et al. 2001). X-ray (Rephaeli, Gruber & 
Persic 1995) and submillimetre (Thomas et al. 2002) de- 
tections have also been reported. The basic properties of 
NGC 5433, obtained from the NASA/IPAC Extragalactic 
Database (NED) unless otherwise noted, are given in Ta- 
ble □ 

Two galaxies in the vicinity of NGC 5433, located 5 
arcmin (CGCG 191-037) and 9 arcmin (KUG 1359+326) 
to the SW, have been identified as possible companions in 
the Uppsala General Catalogue of Galaxies (Nilson 1973). 
More recently, the trio have been listed as WBL 485 by 
White et al. (1999) in their catalogue of nearby poor clus- 
ters of galaxies. On larger scales NGC 5433 has been paired 
with IC 4357, 44 arcmin to the SW and differing in reces- 
sional velocity from NGC 5433 by 40 kms -1 , in the Nearby 
Optical Galaxies (NOG) group assignment of Giuricin et 
al. (2000). Several other faint galaxies, identified using the 
Minnesota Automated Plate Scanner (MAPS) by Cabanela 
(1999) or in the 2 Micron All-Sky Survey (2MASS), are also 
in the field shown in Fig. Two of these have been labelled 
SIS-1 (MAPS-NGP 0.325.0013717) and SIS-2 (MAPS-NGP 
O.271.0297230) for simplicity. There are no galaxies brighter 
than 17th (blue) magnitude in regions of the same size in 
other directions with respect to NGC 5433. Available data 
(from NED) for the galaxies labelled in Fig. are given in 
Table [U 

In this paper we examine the H I distribution and kine- 
matics of NGC 5433, and characterize its environment for 
the first time. Our observing and data reduction techniques 
are given in and details of the radio frequency interfer- 
ence (RFI) excision we performed are in Appendix^] In 
we present our results: an analysis of the HI content and 



1 The VLA is a facility of the National Radio Astronomy Obser- 
vatory (NRAO). The NRAO is a facility of the National Science 
Foundation operated under cooperative agreement by Associated 
Universities, Inc. 



kinematics of NGC 5433 is in and in Eland EHlwc 
characterize the environment of NGC 5433 and the H I in its 
companions. In f^we discuss the prevalence of high-latitude 
HI in edge-on spirals, the characteristics of the NGC 5433 
group and possible origins for the extraplanar emission ob- 
served. A summary of our findings is in JS] 



2 OBSERVATIONS AND DATA REDUCTION 

The data were obtained with the VLA in its C and D con- 
figurations on 2001 July 21, 23 and 1999 May 4, respec- 
tively. NGC 5433 was at the pointing centre of each obser- 
vation, and as such all of the galaxies listed in Table^except 
IC 4357 fell within the instrument's primary beam. During 
each run the galaxy was observed in 25-35 minute intervals, 
separated by 3-5 minute phase calibrator observations. On- 
line Hanning smoothing was applied. The flux calibrator, 
also used as a bandpass calibrator, was observed at both 
the start and end of each observing run. A total of 554 min- 
utes were spent on source in the C configuration, and 142 
minutes in the D configuration. Windy conditions through- 
out the D configuration run caused a number of antennas to 
stow, reducing the amount of data obtained at those times 
by an average of 10 per cent. 

The data from each observing run were calibrated sep- 
arately using the Astronomical Image Processing System 
(aips). Standard flux, phase and bandpass calibration rou- 
tines were applied, and the continuum emission in the re- 
sulting cubes was removed by subtracting linear fits to the 
real and imaginary parts of the visibility data in the line-free 
frequency channels of each bandpass. Following these steps, 
strong RFI was detected in the central 1 or 2 channels of each 
dataset. The culprit was the 1400 MHz band centre (corre- 
sponding to the recessional velocity of NGC 5433) adopted 
for the observations: with this setup, correlated noise from 
the 7th harmonic of the VLA's 200 MHz local oscillator at 
L-band contaminated the shortest east-west baselines in the 
central channels, for which the expected fringe rate is zero 
(Bagri 1996). This RFI was excised by a) removing base- 
lines with east-west projections close to the array centre 
in all channels, and b) clipping contaminated data in the 
visibility domain from infected central channels only. With 
this approach the RFI was excised at the root mean square 
(RMS) noise level a of the infected channels. The main re- 
sults of the RFI removal are that a for the central channel 
of each dataset is ~10 per cent higher than the correspond- 
ing datacube mean, and flux estimates for NGC 5433 from 
the central 2 channels of the C configuration data have an 
associated systematic error on the order of a; all other flux 
measurements and HI maps are unaffected. Details of the 
RFI excision procedure are given in Appendix 1X1 

After calibration, the C configuration datasets were 
combined in the visibility domain and imaged using a va- 
riety of weighting schemes to emphasize HI structures of 
different spatial scales. The C and D configuration data were 
also combined and imaged in a similar manner; we will refer 
to this combination as C+D data. Each data cube was then 
CLEANed (Clark 1980) and corrected for the attenuation of 
the primary beam. Due to the lower sensitivity (relative to 
the C+D maps) and resolution (relative to the C maps) of 
the D configuration data, we show no images of this dataset 
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Figure 1. DSS image of NGC 5433 and its SW environment. Labels have been placed immediately below the galaxy. The more distant 
galaxy IC 4357 is not shown; it is located 44 arcmin SW of NGC 5433, along a line joining NGC 5433 and CGCG 191-037. 



Table 1. Basic Properties of NGC 5433 and Nearby Galaxies 



Galaxy 


Morph. 


a (J2000) 


S (J2000) 


Sep. a 


Mag. 


Diameter 


V b 

v sys 


i c 




Type 


/h m B J 


(° ' ") 


(arcmin) 


(blue) 


(blue; arcmin 2 ) 


(km s ) 


(°) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NGC 5433 


SAB(s)c: d 


14 02 36.05 


32 30 37.8 




14.10 


1.6 x 0.4 


4354 ± 7 


78 


CGCG 191-037 




14 02 20.56 


32 26 53.5 


5.0 


15.6 


0.8 x 0.2 




76 


KUG 1359+326 




14 01 55.87 


32 27 29.6 


9.0 


15.7 


0.4 x 0.2 




60 


SIS-l e 




14 01 49.62 


32 20 34.5 


14.0 


17.45 


0.37 x 0.17 




63 


SIS-2 e 




14 02 38.97 


32 27 50.3 


2.9 


18.00 


0.21 x 0.16 




40 


IC 4357 


S? 


14 00 43.69 


31 53 39.0 


44.0 


14.75 


1.1 x 0.6 


4394 ± 36 


60 



a) Projected angular separation from NGC 5433. 

b) Heliocentric, optical definition (and so throughout). 

c) From Col. 7 assuming a thin disc except for NGC 5433, where an intrinsic axial ratio q = 0.13 is adopted. 

d) Revised from RC3 (de Vaucouleurs et al. 1991) classification of Sdm, by Dale et al. (2000). 

e) Alternates: SIS-1 = MAPS-NGP O.325.0013717, SIS-2 = MAPS-NGP O.271.0297230 



3P 



§ 26 



24 



22 



20 



alone. However, the C and D configuration observations were 
obtained independently and have similar sensitivities (see 
Table |5J, and thus the D data provide a useful consistency 
check on the reality of faint features detected in the C data: 
one expects to find them in both cubes. All maps displayed 
are uncorrected for the primary beam, but calculations are 
performed on the corrected cubes. The observing and map 
parameters for the C, D and C+D data are given in Tabled 

A map of the continuum emission was also made from 
the line-free channels and used to measure continuum fluxes 



in the vicinity of our H I detections (Table . Since the 
resulting images of NGC 5433 do not improve upon those of 
Irwin et al. (1999), we do not reproduce them here. 



3 RESULTS 

In this section we present our findings regarding the H I con- 
tent of NGC 5433 and its environment. In H3.ll we use the 
higher resolution C configuration data to examine in-disc 
and extraplanar HI of NGC 5433, and adopt the higher 
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Table 2. HI Observing and Map Parameters 



Parameter 


C config. 


D config. 


C+D combined 


(1) 


(2) 


(3) 


(4) 




Observing date 


2001 July 21 & 23 


1999 May 4 






On-source observing time (min.) 


554 


142 


696 




Band centre (km s — 1 ) 


4352 


4352 


4352 




Total bandwidth (km s _1 ) 


657 


657 


657 




Channel width (km s _1 ) 


10.1 


21.2 


21.2 




Natural weighting: 










bynth. beam (arcsec) <9 PA ( ) 


19 x 17 6 70 


53 x 50 @ -24 


30 x 29 © 


-54 


RMS map noise a a (mjy beam -1 ) 


0.31 


0.36 


0.21 




Robust weighting 6 : 










Synth, beam (arcsec) @ PA (°) 


16 x 15 & 74 


47 X 45 @ -24 


23 x 21 © 


-58 


RMS map noise a a (mjy beam - 1 ) 


0.32 


0.38 


0.21 




Uniform weighting: 










Synth, beam (arcsec) @ PA (°) 


13 x 13 ® 68 


42 x 39 @ -26 


16 x 15 @ 


-53 


RMS map noise a a (mjy beam -1 ) 


0.37 


0.41 


0.24 





a) a at V = 4352 kms -1 is ~10 per cent higher than this value; see text and Appendix IaI 

b) AIPS robustness parameter is chosen to produce a synthesized beam intermediate to natural 
and uniform weighting. 




Figure 2. Global profiles for NGC 5433 (solid lines), KUG 
1359+326 (long-dashed lines), SIS-1 (short-dashed lines), and 
SIS-2 (dotted lines). Solid error bars represent la random errors. 
For NGC 5433, the dotted part of the error bars represents ad- 
ditional systematic errors on those points from RFI removal; see 
J3]and Appendix IaI for details. 



sensitivity C+D data to characterize the environment of 
NGC 5433 in J£l The C+D data are also used to study the 
HI content of the detected companions in H3.3I Throughout, 
we show images of datacubes with weightings that highlight 
the features we wish to discuss. 



3.1 HI in NGC 5433 

The spatial/spectral resolution and sensitivity of our C con- 
figuration data allow us to examine the morphology and 
kinematics of the HI in NGC 5433 in some detail. We dis- 



cuss its global HI content and parameters in H3.1.1I model 
its H I distribution in H3.1.2I and characterize its extraplanar 
HI extensions in H3.1.3I 



3.1.1 Global H I content 

The global profile of NGC 5433, obtained from the naturally 
weighted C configuration data, is shown in Fig. |2] Channel 
maps are given in Fig. [3] HI is detected from 4586 km s -1 
to 4151 km s _1 , with the north side of the galaxy receding 
and the south side advancing. Both figures show that HI 
emission is stronger on the southern side of the galaxy. This 
asymmetry is also seen in the single-dish global profiles of 
Mirabel & Sanders (1988) and Schneider et al. (1990). Note 
that the "dip" in the global profile near 4350 km s~ may 
not be real, given the larger (dotted) error bars in these 
channels (see f|5|and Appendix IaI. 

Fig. [I] shows the total intensity (zeroth moment) and 
intensity- weighted velocity field (first moment) of NGC 5433 
from uniformly weighted C configuration data. A mild warp 
is present, most easily seen toward the ends of the opti- 
cal disc. This warp is complex, showing curvature toward 
the west (east) at the northern (southern) edges of the op- 
tical disc and in the opposite direction (east on the north 
side, west on the south side) farther out. 3 features extend- 
ing from the disc, labelled Fl through F3, are also visible. 
From the mean column densities and sizes of the features 
in Fig. we estimate that an H I mass Mhi ~ 10 s Mq is 
associated with each. The velocity field (Fig. 0|d) is well or- 
dered, but the kinematic major and minor axes are not quite 
perpendicular. This may indicate the presence of a bar or 
other oval distortion in the disc, and is perhaps expected 
given NGC 5433 's morphological classification as a barred 
galaxy (Dale et al. 2000). The extensions Fl - F3 show no 
kinematically distinct signatures in the first moment map. 

A position-velocity plot for a 3 arcsec-wide slice from 
the uniformly weighted data along the major axis of 
NGC 5433 is given in Fig. |^| It reveals that the disc ex- 
hibits largely solid-body rotation out to 20 arcsec (4 kpc, 
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Figure 3. Robustly weighted C configuration channel maps for NGC 5433. Contours are at 0.31 X (-3 (dashed), -2 (dashed), 2 (2cr), 3, 
5, 7, 10, 14, 18, 22, 26, 30, 32, 34, 36, 37, 37.5) mjybeam - 1 . The channel velocity in kms - 1 is indicated at the top of each frame, and 
the synthesized beam is given in the lower left corner of the first frame. 
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Figure 3. Continued. 



after correction for beam smearing) from the centre. The 
rotation curve of the galaxy flattens out only in the H I-rich, 
southern side of the disc. From this figure as well as Fig.^K, 
it is clear that the northern HI does not extend as far out 
as the southern gas. 

In Table [3] we present global properties of NGC 5433 
obtained from the C array maps. The systemic velocity V sys , 
measured line width at 50 per cent of the peak W sin i, 
and integrated flux JSdV are in rough agreement with 
previously published values from single-dish measurements 
(Mirabel & Sanders 1988, Schneider et al. 1990). This sug- 
gests that we have recovered all of the H I emission and we 



thus make no short-spacing corrections to JS dV or in our 
galaxy models (see 33.1.21 . The HI mass is computed as- 
suming an optically thin disc as in Shostak (1978), where D 
is the distance to NGC 5433 in Mpc: 



M H i = 2.36xl0 5 D 2 



SdV Mq. 



(1) 



The total dynamical mass is computed assuming a spherical 
dark matter distribution: 

M T = (6.78 x 10 4 )6 r DV 2 (r) M . (2) 
In Equation 6 r is the disc radius in arcmin and V(r) is 
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Table 3. Global Properties of the Galaxies 



Gcilfixy 

(i) 


T/ a 
v sys 

(kms- 1 ) 
(2) 


Wo sin % 
(km s _1 ) 
(3) 


JbdV 
(Jykms 1 ) 
(4) 


(x 10 9 M Q ) 
(5) 


IVl'J' 

( x 10 11 M ) 
(6) 


J20 
(mjy) 
(7) 


NGC 5433 d 
KUG 1359+326 
SIS-1 
SIS-2 


4346 ± 7 
4282 ± 29 
4321 ± 16 
4431 ± 19 


381 ± 7 
123 ± 29 
163 ± 16 
91 ± 19 


5.89 ±o;ols 
1.01 ± 0.11 
0.66 ± 0.07 
0.12 ± 0.04 


5-89 ± il r s 

1.01 ± 0.11 
0.66 ± 0.07 
0.12 ± 0.04 


1.64 ± 0.14 
0.10 ± 0.03 
0.16 ± 0.02 
0.05 ± 0.01 


62.8 ± 3.2 
2.0 ± 0.4 

< 0.59 

< 0.36 



a) Midpoint of the global profile (Fig. |5J at 50% of the peak. 

b) Velocity width at 50% of the global profile peak (Fig.[2J. 

c) Continuum flux densities. Upper limits (3cr) assume an unresolved point source. 

d) For NGC 5433, quantities derived from total channel fluxes have both random and systematic er- 
rors (denoted r and s, respectively), the latter a result of RFI excision. See J5]and Appendix^Jfor details. 



the inclination-corrected rotation velocity. We make no cor- 
rection for turbulence, as it is expected to be much smaller 
than the quoted uncertainties (Tully & Fouque 1985). For 
NGC 5433, we estimate Mr for the approaching and reced- 
ing sides of the galaxy separately adopting as 9 r and V(r) 
the lowest contour in Fig.0 and average the result to obtain 
the value in Table |3] 

In light of the recent morphological re-classification of 
NGC 5433 to SAB(s)c: by Dale et al. (2000) from the origi- 
nal RC3 classification of Sdm, it is interesting to compare its 
global properties with those of a sample of normal field spi- 
rals of various morphological types. For this exercise we use 
the volume-limited RC3-LSc sample of Roberts & Haynes 
(1994) as the comparison sample, and compute the physi- 
cal parameters of NGC 5433 using the single-dish data of 
Schneider et al. (1990) for the HI parameters (the results 
are the same if our synthesis data are used). The results in 
Table [I] show that, although there is overlap in properties 
between the morphological types, the physical parameters 
of NGC 5433 are indeed more typical of the newer Sc clas- 
sification than for Sd or Sm. 



3.1.2 Kinematic modelling 

Our H I data for NGC 5433 contain intensity measurements 
at many positions and velocities above the noise, all of which 
provide important information about the galaxy's global H I 
distribution and dynamics. Thus, rather than modelling the 
moment maps alone, we have modelled the HI spectra at 
all positions in the cube. This is done by parameterizing the 
velocity and density distributions in the galaxy, specifying a 
geometry, producing model HI spectra, and then searching 
through parameter space for the best fit. The benefits of this 
method are that all data in the cube are used, the volume 
density distributions (rather than just the column densities) 
and the system geometry can be constrained, and it is pos- 
sible to disentangle the effect of a thick disc from the effect 
of inclination on the projected distribution perpendicular 
to the major axis. In addition, examination of the residual 
cube can provide more insight into the peculiarities of the 
H I distribution in the galaxy. 

The routine used (called CUBIT and freely avail- 
able from http : / /www . astro . queensu . ca/ "irwin) adopts a 
Brandt rotation curve and can also apply a velocity dis- 
persion if it is warranted. The density distribution in the 
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Figure 5. 3 arcsec-wide position-velocity slice of the uniformly 
weighted data shown in Fig. [3] along the major axis of NGC 5433. 
Contours are at (-1, 1, 2, 3, 4, 5, 6, 7, 8, 9) mjybeam -1 , and 
greyscale runs from to 9.08 mjybeam -1 . Coordinates are rela- 
tive to optical centre and recessional velocity of NGC 5433 (Ta- 
bic □ . 



disc can be either constant, decreasing exponentially, gaus- 
sian, or a gaussian ring with a centre at a specified radius 
and different scale lengths on either side of this position. 
The latter distribution is useful in modelling galaxies with a 
dearth of H I at small galactocentric radii. The most recent 
version of the routine also allows for two radial distributions 
to be superimposed. The vertical density distribution can be 
either constant, exponential, or gaussian. Details of this pro- 
cess can be found in Irwin & Chaves (2003) and references 
therein, or at http://www.astro.queensu.ca/ "irwin. 

For NGC 5433, we modelled the C configuration data at 
three different weightings/resolutions (natural, robust and 
uniform; see Table |5J to ensure that the results are inde- 
pendent of the imaging parameters. This covers a range of 
spatial resolutions from 13 - 19 arcsec (4.1 - 5.7 kpc). The 
quoted errors on each best-fitting parameter in the model 
reflect the variation in results over the variously weighted 



© 2004 RAS, MNRAS 000.HHTS1 



8 Spekkens, Irwin and Saikia 




II 




4520 -( 
- 4400 








4300 ^ — x\ 


- 


4200 




(b) 




ii 





38 37 36 35 34 
RIGHT ASCENSION (J2000) 



RIGHT ASCENSION (J2000) 



Figure 4. (a) Total intensity map of NGC 5433 (contours; from C array uniformly weighted data) overlaid on an optical image (greyscale; 
DSS). Contours are at 63.5 X (2, 3, 4, 5, 6, 7, 8, 9, 10) Jy m s _1 beam _1 . These values would be multiplied by 6.70 X 10 18 to convert to 
column density, Njji, in units of cm -2 , (b) First moment map of NGC 5433, from the C data. Contours are in 20 kms -1 increments. 



Table 4. NGC 5433 and Typical Properties of Late-Type Spirals 



Galaxy 
(1) 



log(R lin ) b 
(2) 



log(Lg) 

(3) 



log(Mr) 
(4) 



logOflr) c 
(5) 



log( M H i) 
(6) 



\og(M HI /L B ) 
(7) 



NGC 5433° 


1.18 ± 0.05 


10 


.30 


11.17 


± 0.05 


0.84 ± 0.09 


9.69 


± 0.05 


-0.61 ± 0.05 


Sc d 


1.13- 


+ 1.40 


10.2- 


-+10.9 


10.8- 


+ 11.3 


0.80- 


+ 1.27 


9.8- 


^10.4 


-0.3- 


+ -0.7 


Sd d 


0.85- 


+ 1.22 


9.8- 


+ 10.3 


10.2- 


-+10.7 


0.83- 


+ 1.22 


9.4- 


^10.0 


0.0- 


► -0.5 


Sm d 


0.88- 


+ 1.16 


9.4- 


+ 10.2 


9.9- 


+ 10.4 


0.82- 


+ 1.28 


9.3 


^9.8 


0.1- 


► -0.3 



a) Parameters as in Roberts & Haynes (1994), using single-dish data from Schneider et al. (1990). 

b) Linear radius (kpc) from RC3. 

c) HI surface density in Mq pc — 2 , defined as cfhi = Mhi /^Rf in - 

d) Range of values (25th to 75th percentile) from RC3-LSc sample of Roberts &: Haynes (1994), adjusted 
to Ho = 70 kms -1 Mpc -1 . 



cubes. Although there are up to 15 free parameters in the 
model, there are over 250 independent data points above 3a 
in the robustly weighted cube; thus the parameters are well 
constrained by the data. 

The best-fitting model values for NGC 5433 are given in 
Tableland moment maps made from the residuals (data mi- 
nus model) are displayed in Fig.|S|for the robustly weighted 
cube. Model 1 gives the best results for the case in which 
a single radial density distribution is used, and Model 2 is 
the best fit for the sum of two radial density distributions. 
For both models the RMS noise of the residual cube is less 
than 3cr of the input data, suggesting that the global mor- 
phology and kinematics of NGC 5433 are well-represented 
by the models. The forms of these distributions are given in 
the notes to Table |K| 

The first 5 parameters of Table [S] specify the kinematic 
centre of the galaxy and the disc geometry. V sy s agrees with 



the value from the global profile (Fig. |2J but the kinematic 
centre is 3.3 arcsec south of the optical centre, likely due 
to the observed HI asymmetry with more gas weighted 
toward the south (Fig. QJ. The remaining parameters de- 
scribe the shape of the rotation curve and the density dis- 
tributions within and perpendicular to the midplane. In the 
midplane, the best-fitting distribution is a ring centred at 
Ro = 12.6 kpc. The RMS of the residual cubes, in the last 
row of Table |S] shows that the addition of the second radial 
density distribution peaking at the centre with a character- 
istic radius of D c — 3.5 kpc (Model 2) improves the fit only 
marginally. Our models therefore do not provide conclusive 
evidence for a central HI depression in NGC 5433. For both 
Model 1 and Model 2, the vertical exponential scale height 
(H e ^ 1.2 kpc) is an upper limit because the result is com- 
parable to the pixel size adopted for the cube. 

The moment maps of the residual emission (Fig. |SJ 
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Table 5. Best-Fitting Kinematic and Density Parameters 



Parameter 


Model 1 


Model 2 




(1 Radial Dist.) 


(2 Radial Dist.) 


(1) 


(2) 


(3) 


ni a ( 12000") (h m s") 


1 4 02 SB 00 + 05 


1 4 02 35 99 + 06 


8 a (J2000) (° ' ") 


32 30 34.5 ± 0.5 


32 30 34.5 ± 0.5 


PA 6 (°) 


3.5 ± 0.4 


3.1 ± 0.4 


i (°) 


81.0 ± 0.2 


83 ± 2 


Vc,c (km R ^ 
v sys \" ' ' ' a f 


4356 ± 3 


4355 ± 3 


V mn ^ c f+m s ^ 


198.8 ± 0.2 


198 ± 1 


R „, „ _ c (knrl 
L **vnax V P W 


6 + 1 


5.3 ± 0.3 




2.0 ± 0.6 


1.5 ± 0.4 


no d (cm -3 ) 


0.139 ± 0.007 


0.13 ± 0.01 


R d (kpc) 


12.6 ± 0.3 


12.6 ± 0.3 


D D d (kpc) 


2.18 ± 0.09 


2.24 ± 0.03 


Di d (kpc) 


10.1 ± 0.6 


10.1 ± 0.6 


D c d (kpc) 




3.5 ± 0.3 


n c ™ (cm -3 ) 




0.11 ± 0.01 


H e e (kpc) 


<C 1.2 ± 0.1 


^ 1.2 ± 0.3 


RMS^ (mjy beam" 1 ) 


0.49 


0.45 



a) Kinematic centre of the H I distribution. 

b) PA of the receding major axis. 

c) Rotation curve parameters, defined by the Brandt curve: 
V(R) = Vmax{R/Rma^)[l/3 + 2/3 (R/R max )™] <" 3 / 2m ) 
where R is the radial distance in the plane of the galaxy and 
V(R) is the rotational velocity at that radius. 

d) The radial volume density distribution parameters are given 
by n(R) = ni(R) + n2(R), where 

ni {R) = n exp[-{R - Ro) 2 /D 2 )] for R> R , 
ni (R) = n exp[-(R () - R) 2 /Di 2 )} for R sC R , 
n 2 {R) = n c exp[~(R 2 /D 2 )}. 

e) Perpendicular exponential scale height. The density at any 
point in the galaxy is n(R, z) = n(R) exp(— z / H e ) , where z is 
the perpendicular distance from the plane of the galaxy. 

/) RMS of the residual (model-data) cube after applying a 2a 
cutoff, where a is the RMS noise level of the original input data. 



clearly show the HI asymmetry mentioned previously, with 
excess emission on the south side of the galaxy. The remain- 
ing features in the residual emission appear to be related to 
the extensions, discussed further in the next section. Note 
that these residuals can not be eliminated by varying the 
model parameters, such as the input position angle (PA). 



3.1.3 High-latitude HI in NGC 5433 

The residual maps in Fig. EJreveal departures from the best- 
fitting models of NGC 5433 at high latitudes, coincident 
with the features Fl - F3 in the total intensity maps (Figs. 2] 
and|SJi). We examine these extensions in the C configuration 
data, which is most suited to the analysis given its spatial 
and spectral resolution. Since NGC 5433 is not truly edge- 
on (i — 78°, Table0 and since many features of interest are 
comparable in size to the map resolution, the linear dimen- 
sions we infer from projected angular sizes in the map are 
only approximations to the true physical scales involved. 

In Fig. [TJl, we show a total intensity H I map of the ro- 
bustly weighted datacube for NGC 5433, enhanced to illus- 
trate low intensity extensions, with the high resolution radio 
continuum emission contours of Irwin et al. (2000) superim- 
posed. Extensions from the mid-plane are clearly present in 




14 02 39 38 37 36 35 34 33 39 38 37 36 35 34 33 
RIGHT ASCENSION (J2000) 

Figure 6. Moment maps of robustly weighted NGC 5433 data 
and residuals of corresponding models. Total intensity contours 
may be converted to column densities Nhi (in units of cm -2 ) by 
multiplying levels by 4.5 X 10 18 . (a) Integrated intensity (zeroth 
moment map) of the data. Contours are at 10 X (2.5, 5, 10, 20, 
30, 40, 50, 60, 70, 80) Jy ms -1 beam -1 . The features Fl - F3 
identified in Fig. |1] are labelled, (b) First moment map of the 
data, (c) Zeroth moment map made from the residual (data - 
model) cube for Model 1. Contours are at 10 X (2.5, 5, 10, 20) 
Jyms _1 beam _1 . (d) First moment map of the data in 6b minus 
the first moment map made from Model 1. The greyscale ranges 
from -50 (light) to +50 (dark) km s -1 and contours are at -20 
km s _1 (dashed) and +20 km s — 1 (dotted), (e) Same as in 6c, 
but for Model 2. (f) Same as in 6d, but for Model 2. 



both ISM tracers. As has been observed in other galaxies 
(e.g. Collins et al. 2000), when extensions exist, they can 
be seen in a variety of components. The low intensity emis- 
sion shows that the HI features extend as far as 25 arcsec 
(~5 kpc) in projection from the midplane. This is greater 
than the disc thickness 1.2 kpc, t!3.1.21 but appears to be 
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typical of galaxies showing HI extensions (e.g. Lee & Irwin 
1997). Feature F3, which appears as a broad single feature 
in Fig. is revealed to consist of at least two components, 
one to the north and one to the south of the nucleus. The 
continuum image shows two similar features on the west 
side, despite its drastically different resolution of 1.5 arcsec 
compared to that of the H I data. Given the correspondence 
of the western extensions in the maps, it may be that each 
continuum feature is related to the nearest H I extension and 
possibly to each other as well. 

It is not always clear whether HI extensions identi- 
fied in total intensity maps represent coherent structures, or 
whether they consist of unrelated features along the same 
line of sight. To distinguish between these possibilities, we 
have searched through the data both in position-position 
(P-P) space as well as in position-velocity (P-V) space to 
look for coherent structures related to the features Fl, F2, 
and F3. 

In the case of Fl (14 h 02 m 38 s , 32° 31' 15"), 2 - 5a exten- 
sions from the underlying disc are seen at 4522 and 4490 - 
4469 kms" 1 (Fig.0. The features from 4490 - 4469 kms' 1 
are also evident in the C+D data, and given the direction 
of rotation of NGC 5433, their NE orientation is consistent 
with the slight redshift of the data relative to the best-fitting 
disc models in this vicinity (Figs. EJi and EH). If they repre- 
sent a single, coherent structure, then using the fluxes mea- 
sured for these extensions over the 4490 - 4469 kms" 1 ve- 
locity range in Equation Q we compute Mhi ~ 7x 10 7 Mq. 
This is broadly consistent with Mhi estimated for Fl from 
the total intensity distribution in Fig. [IJl 

At the position of F2 (14 h 02 m 37 s 5, 32° 30' 40") there 
are many features in the channel maps, such as 2 - 3cr peaks 
from 4447 - 4426 kms" 1 and a series of extensions from 
4331 - 4299 kms" 1 . However, there is no evidence that these 
features form a coherent structure. This is illustrated in the 
P-V slice shown in Fig. |7j3: near F2 (corresponding to a 
vertical strip to the left of the main emission), there are a 
number of small, likely unrelated extensions at a variety of 
velocities which collectively contribute to F2. 

F3 is located at 14 h 02 m 34 s , 32° 30' 35" (approximate 
midpoint). We see 2 - 3<r extensions in Fig.|3]at that location 
at 4437, 4405, 4394 and 4352 kms" 1 , and 2 - 5a ones from 
4256 - 4235 kms" 1 . The latter are also present in the C+D 
data, and are consistent with the residual maps in Figs |HJi 
andHJ. Considering features in this velocity range to repre- 
sent a coherent structure, we compute Mhi ~ 1.1 x 10 8 Mq 
from Equation (1), again consistent with the mass computed 
for F3 in Fig. EJ,. 

The P-V slice through the top part of F3 (Fig. 0J>, 
right side of main emission) further illustrates that the 
emission does appear to form a coherent structure which 
is roughly circular, with the centre at a~ 14 h 02 m 35 s , V ~ 
4250 kms" 1 : this is a signature of an expanding shell or part 
thereof. The diameter of the shell-like feature is 50 km s" 1 in 
velocity and 25 arcsec (~3 kpc) in right ascension. An esti- 
mate of the H I mass from this P-V slice only is ~ 3 x 10 7 Mq . 
Under the assumption that the shell is expanding, we thus 
compute an expansion kinetic energy ~ 2 x 10 53 ergs, clearly 
a lower limit for that associated with the feature as not all 
the mass is represented in this slice. A P-V slice through 
the lower part of F3 (Fig.[7p) shows more emission at a sim- 
ilar velocity, and is typical of what is seen in other slices 



through this feature. If all of the mass computed for this co- 
herent part of F3 (1.1 x 10 s Mq) is involved in the expansion, 
then the associated kinetic energy of the HI is ~ 8 x 10 53 
ergs. 

If the feature in Fig. [TJi is generated from a shperically 
symmetric, point injection of energy such as a supernova, 
then the required input energy can be computed using the 
results of numerical simulations by Chevalier (1974) for an 
expanding shell that is now in the radiating phase of its 
evolution: 



C O w 1(1*3 1.12 D 3.12 rr 1.4 

E = 5.3 x 10 m Rsh Vsh 



ergs, 



(3) 



where V s h (km s" 1 ) and Rsh (pc) are the expansion velocity 
and radius of the shell, and ri\ (cm" 3 ) is the ambient den- 
sity at the location of generation. Using V s h and R s h of the 
structure in Fig|7J) and estimating ni at its location from 
the models in g3X2] then from Equation @£~2x 10 54 
ergs. This estimate decreases only by a factor of 2 if, instead, 
a continuous deposition of energy is assumed (Weaver et al. 
1977). 

One alternative is to interpret the apparently circular 
structure of Fig. [7Ji as disconnected components, pushed 
or pulled out of the midplane and collectively contribut- 
ing to F3. However, even if the clouds aren't related one 
still needs to place ^ M c i ou d ~ 10 8 Mq an average distance 
z ~ 2.5 kpc above the plane. If zo is the vertical scale height 
of the stellar mass distribution and po is the stellar mass den- 
sity at the midplane, then the potential energy of a "cloud" 
in equilibrium at a distance z above a thin, axisymmetric 
disc is 



$(z) =2xl0 54 (- 



Mdoud 



) (-=*-) ( 



Po 



2 x 10 7 M Q ' v 700pc' v 0.185M Q pc- 



r) ln [ cosh (-)] er e 



(c.f. Binney & Tremaine 1987, Lee et al. 2001), where we 
normalize zq to Galactic and po to solar neighbourhood 
values, respectively. Adopting these values for zq and po 
we find $ ~ 3 x 10 55 ergs for £ M cioud = 10 8 M and 
z — 2.5 kpc, larger than the injection energy computed 
above. While these energy estimates are order-of-magnitude 
approximations (c.f. Lee & Irwin 1997), it is clear that the 
energy required to produce F3 is very high, and comparable 
to inferred values for the largest supershells in the Galaxy 
(McClure-Griffiths et al. 2002) and those in other edge-on 
spirals (e.g. Chaves & Irwin 2001). 



3.2 The environment of NGC 5433 

To characterize the environment of NGC 5433, we imaged 
the C+D data with a variety of weighting schemes to em- 
phasize structures on different spatial scales, and system- 
atically searched the resulting HI cubes. In all, there are 
20 galaxies within a 16 arcmin radius (the primary beam 
HWHM) of NGC 5433 with MAPS or 2MASS identifiers. 
Of these, we detect HI in 3 galaxies other than NGC 5433 
itself: KUG 1359+326, SIS-1 and SIS-2 (see Fig.[TJ. We note 
again here that IC 4357, a confirmed companion 44 arcmin 
SW of NGC 5433 (Giuricin et al. 2000), is not within the 
primary beam of our observations; we do not search for emis- 
sion outside this region. We find no other candidates with 
peak fluxes over 4a spanning more than 2 consecutive chan- 
nels. Adopting these parameters as detection criteria, our 
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Figure 7. (a) High resolution continuum map of NGC 5433 (contours) taken from Irwin et al. (2000), superimposed on an HI zeroth 
moment map (greyscale, same as Fig. 6a) with the greyscale range from 50 to 250 Jy m s -1 beam -1 . Horizontal lines show locations of 
P-V slices in 7b and 7c. (b) P-V cut taken along the top horizontal line in 7a. Contours are at 0.6, 0.75, 1.2, 2.0, and 4.0 mjy beam 
and the greyscale is from 0.3 to 3.9 mjy beam -1 , (c) Same as 7b, but for a P-V cut along the bottom horizontal line in 7a. 



maps are sensitive to Mhi = 4 x 10 7 Mq at the distance of 
NGC 5433. 

The detection of H I in KUG 1359+326 (projected sep- 
aration = 200 kpc), confirms that this galaxy is indeed 
a physical companion to NGC 5433. We do not detect 
CGCG 191-037, despite its proximity to NGC 5433 (pro- 
jected separation = 102 kpc), diameter intermediate to 
NGC 5433 and the companion KUG 1359+326, and blue 
magnitude comparable to the latter (Table 1). Its near- 
infrared brightness is also intermediate to that of these two 
group members, as it is an average of -~1.25 magnitudes 
fainter than NGC 5433 and ~1.27 magnitudes brighter than 



KUG 1359+326 across the 2MASS J, H and K bands. Per- 
haps this galaxy's recessional velocity places it outside our 
657 kms -1 bandpass centred on NGC 5433. In this sce- 
nario CGCG 191-037 is at least 4.5 Mpc from NGC 5433, if 
their relative velocity is dominated by Hubble flow. While 
the size and magnitudes of CGCG 191-037 are compara- 
ble to those of KUG 1359+326, its magnitude relative to 
the dominant group member NGC 5433 suggest that it is a 
background (rather than a foreground) object. Alternately, 
CGCG 191-037 may have a substantial peculiar velocity rel- 
ative to NGC 5433, or may be devoid of HI altogether. 
We also searched for continuum emission in CGCG 191- 
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Figure 8. Naturally weighted C+D channel maps for KUG 1359+326 (in the upper half of the frames) and SIS-1 (in the lower half). 
Contours are at 0.21 X (-3 (dashed), -2 (dashed), 2 {2a), 3, 4.5, 6, 8, 10, 12, 14, 16) mjybeam - 1 . The channel velocity in kms - 1 is 
indicated in the upper right corner of each frame, and the synthesized beam is given in the lower left corner of the first frame. 
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Figure 9. Naturally weighted C+D channel maps for SIS-2. Con- 
tours arc at 0.21 X (-3 (dashed), -2 (dashed), 2 (2a), 3, 4, 5, 6, 
7) mjybeam - 1 . The channel velocity in kms" 1 is indicated in 
the upper right corner of each frame, and the synthesized beam 
is given in the lower left corner of the first frame. 



037, and found none. If our non-detection of CGCG 191-037 
in HI stems from this galaxy being closer or farther from 
NGC 5433 than its optical dimensions suggest, the contin- 
uum non-detection supports the latter possibility. 

Particularly interesting is the detection of two new com- 
panions, SIS-1 and SIS-2. SIS-1 is a 17.5th magnitude galaxy 
with a discernible disc, and its projected separation from 
NGC 5433 is large at 289 kpc. SIS-2 is an 18th magnitude 
galaxy just south of NGC 5433, separated from its centre by 
2.9 arcmin (55 kpc). It is clear that NGC 5433 is the domi- 
nant member of a group of at least 4 galaxies (5 if IC 4357 
is included). At the sensitivity of our observations, all the 
H I companions within 300 kpc of NGC 5433 are to the SW 
of the latter. 



3.3 HI in the companions 

Although the companions of NGC 5433 are detected in both 
our C and C+D maps, we analyse of their H I distributions 
using the higher sensitivity C+D data. Naturally weighted 
C+D channel maps for KUG 1359+326, SIS-1 and SIS-2 
are shown in Figs. |H1 and and the global profiles derived 
from these maps are found in Fig. [5] Moment maps for each 
galaxy, again obtained from the naturally weighted C+D 
data, are given in Fig. 1101 The global HI properties are 
listed in Table [3] In computing Mhi and Mt from Equa- 
tions Q and we assume that the companions are at 
the same distance as NGC 5433. For Mt, we approximate 
8 r = 1.5 x a/2 where a is the optical major axis diameter of 
each system (Table 0. Given the resolution of these obser- 
vations and kinematic distortions observed (see below), Mt 
only provides an indication of the true dynamical masses of 
the companions. 

We discuss the H I properties of each galaxy below. Be- 
cause of the low spatial and spectral resolution of the detec- 
tions, we do not attempt to model their HI distributions. 



3.3.1 KUG 1359+326 

Our strongest detection among the companions, 
KUG 1359+326 shows H I emission from 4373 - 4204 kms" 1 
(Fig. |8J, and is resolved in most of these channels. The NW 
side of the galaxy is receding, and the SE side advancing. 
The most interesting feature of the maps is an extension 
apparently emanating from the SE part of the disc, primar- 
ily from 4331 - 4288 kms" 1 . It extends ~ 2 arcmin from 
the optical centre of the galaxy, and farther if the low-level 
emission at 4267 and 4182 kms" 1 is also associated with 
the extension. Despite the low S/N of this feature, we detect 
it in both the C and D data independently: this strongly 
suggests that it is indeed real, rather than a chance thermal 
fluctuation in one of the datasets. The extension is also 
seen at a low level in the total intensity map in Fig. 110b . 
At the distance NGC 5433 (and after correcting for beam 
smearing), emission is found over 30 kpc away from the 
disc of KUG 1359+326, over 4 times the optical major axis 
diameter of the galaxy. There is no coherent velocity field 
associated with the extension in our data, and as such it is 
not shown in Fig. 110b . The origin of this extended feature 
is thus unclear from these data. 

The first moment map of KUG 1359+326 (Fig. ITDb ) 
indicates that the velocity field of this galaxy is highly dis- 
torted. Most of the gas from 4270 - 4300 kms -1 appears to 
lie along the optical major axis of the system, although the 
shape of the isovelocity contours indicates the presence of 
kinematic distortions along the minor axis at these veloci- 
ties. The emission at the NE edge of the galaxy, coincident 
with a low-level extension in the same direction in Fig. 110b . 
shows the most distorted velocity contours, which appear 
to be 'pushed' northwards relative to those expected for an 
inclined spiral. 



3.3.2 SIS-1 

The channel maps for SIS-1 (Fig. |SJ show emission from 
4394 - 4246 km s -1 , with the NW side of the galaxy advanc- 
ing and the SE side receding. The low-level emission in the 
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channel maps is barely resolved by our synthesized beam, 
and shows a number of distortions. Two of these features, 
extending to the NE and NW of the emission centroid pri- 
marily at 4352 km s _1 , are also evident in the total intensity 
maps m Fig. [TUt- The NW feature appears as an extension 
from the tip of the major axis of the galaxy. The higher 
level contours in the zeroth moment map are slightly elon- 
gated along the optical major axis of SIS-1, their centroids 
roughly coincident with the optical peak. The velocity field 
in Fig. I10H shows ordered rotation in the SE part of the disc, 
with the kinematic major and minor axes roughly aligned 
with the optical ones. The NW side appears to be kinemati- 
cally distinct from the remainder of the disc, and associated 
with the NW feature of the total intensity map. 

3.3.3 SIS-2 

SIS-2 is the weakest source that we detect in the vicinity of 
NGC 5433, extending from 4479 - 4373 km s _1 in Fig. El It 
is none the less a definite detection in our data, with peak 
emission greater than 6cr over four consecutive channels in 
the C+D maps as well as independent detections in the C 
and D datasets. The north side of the galaxy is advancing 
and the south side receding, and thus its rotation is retro- 
grade relative to that of NGC 5433. We do not resolve the 
emission in most of the channels. As such, the total intensity 
map in Fig. 110b is mostly featureless except for an offset in 
the H I emission peak from the optical centre by about 7 arc- 
sec, or 1/6 of a beam. The velocity field of SIS-2 (Fig. EI) 
shows largely ordered, solid-body rotation throughout the 
disc, likely a result of the poorly resolved emission regions 
in each channel. 



4 DISCUSSION 

Our observations probe the HI morphology and kinemat- 
ics of NGC 5433 and its environment with unprecedented 
detail. In light of our results we discuss the prevalence of 
high- latitude HI in edge-on spirals in H4.ll characterize the 
NGC 5433 system as an interacting group of galaxies in H4.2I 
and underscore the possible importance of environmental ef- 
fects in driving the disc-halo phenomenon in H4.3I 

4.1 High-latitude HI in NGC 5433 and other 
spirals 

We identify 3 extraplanar features in NGC 5433 by exam- 
ining moments of the HI distribution and best-fitting kine- 
matic models. 2 of these (Fl and F3) appear to be associ- 
ated with coherent features in P-P and P-V space (Figs. El 
UJ; we find a complete P-V loop in the upper part of fea- 
ture F3 (Fig. 13)), which may represent an expanding shell. 
For a point injection of energy, the implied input energy 
for F3 is E ~ 2 x 10 54 ergs (Equation J3J). This is nearly 
an order of magnitude larger than the corresponding value 
E = 4.6 ± 2.0 x 10 53 ergs for the most energetic new su- 
pershell (GSH 292-01+55) found in the Southern Galactic 
Plane Survey (McClure-Griffiths et al. 2002) , but lower than 
E inferred for supershells found in other spirals, such as fea- 
ture Fl in NGC 2613 (E > 1.4 x 10 55 ergs, Chaves & Irwin 
2001) and F2 in NGC 5775 (£ ~ 2.1 X 10 55 ergs, Lee et al. 



2001). Given the order-of-magnitude uncertainties on these 
estimates (Lee & Irwin 1997), it is clear that the implied 
energetics for the feature F3 in NGC 5433 are comparable 
to those of the largest structures found in the Milky Way as 
well as features in other spirals. 

While the sample remains small, discrete HI features 
and /or thick HI discs have now been observed in galax- 
ies with a wide range of star formation rates (e.g. higher: 
NGC 891, Swaters, Sancisi & van der Hulst 1997; NGC 2403, 
Schaap et al. 2000; low: UGC 7321, Matthews & Wood 2003) 
and in a variety of environments (e.g. isolated: NGC 3044, 
Lee & Irwin 1997; NGC 3556, King & Irwin 1997; compan- 
ions: NGC 5433; interacting: NGC 5775; Irwin 1994). Re- 
gardless of its origin, coherent extraplanar H I emission may 
thus be common in spiral galaxies, as is the case for high- 
latitude cosmic-ray haloes (Irwin et al. 1999). The existence 
of these features implies significant energy transfers from 
the disc into the halo, by mechanisms that remain poorly 
understood. We return to plausible formation scenarios for 
features in NGC 5433 in AO 

4.2 An active environment 

NGC 5433 is in a richer environment than previ- 
ously thought, having at least 3 physical companions 
(KUG 1359+326, SIS-1, SIS-2) rather than the 2 in- 
ferred through proximity and size alone (KUG 1359+326, 
CGCG 191-037). We suspect that other faint MAPS galax- 
ies in the field might also be detectable in HI with obser- 
vations of greater sensitivity. Since we do not detect either 
HI or 20 cm continuum emission from CGCG 191-037, its 
relationship to the detected galaxies remains unclear; HI 
observations with a wider bandpass would resolve this issue. 
If the more distant NOG companion IC 4357 is included, 
the group size increases to 5; HI observations of IC 4357 of 
similar sensitivity to those presented here may also reveal 
physical companions, making the group richer still. We note 
that the X-ray detection reported for NGC 5433 (Rephaeli 
et al. 1995) may represent hot intergalactic gas, sometimes 
found in poor groups (Fukazawa et al. 2002). 

It is interesting that all confirmed group members differ 
in velocity from NGC 5433 by no more than 70 km s _1 
(Table El- SIS-1, for example, at a projected separation of 
248 kpc, is separated in velocity from NGC 5433 by only 
25 kms -1 . Thus the projected separation of the galaxies 
may be close to their true separation. In addition, within 
20 arcmin of NGC 5433 there are over 3 times as many 
2MASS galaxies (12 vs. 4) and 2 times as many optically 
detected galaxies (14 vs. 7) in the NE+SW quadrants than 
in the NW+SE quadrants. We thus speculate that the NE - 
SW orientation of the galaxies in this study may represent 
the largest systems in a "true" filament at least 770 kpc 
in size (the separation between NGC 5433 and IC 4357), 
with little variation in depth. Thus the group geometry is in 
qualitative agreement with galaxy structures on larger scales 
found in surveys (e.g. Cross et al. 2001) and simulations of 
hierarchical structure formation via gravitational collapse 
(e.g. Bertschinger 1998). 

There is evidence that interactions are taking place in 
the group, either among the detected galaxies or with a hot 
intergalactic medium. NGC 5433, itself, shows some optical 
curvature (Fig. and has an asymmetric H I distribution 
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Figure 10. Total intensity HI maps of detected galaxies superimposed on optical DSS images (left hand panels) and their corresponding 
first moment maps (right hand panels). Zeroth moment map contour values would be multiplied by 1.31 X 10 18 to convert to column den- 
sity, Njji, in units of cm -2 . All first moment map contours are separated by 15 km s _1 unless otherwise labelled, (a) Total intensity map 
of KUG 1359+326. Contours are at 39.4 X (0.5, 1, 2, 3, 4, 5, 6, 7, 8, 9, 9.5) Jy ms- 1 beam _1 . (b) First moment map of KUG 1359+326. 
(c) Total intensity map of SIS-1. Contours are at 24.9 X (1, 2, 3, 4, 5, 6, 7, 8, 9, 10) Jy m s-'beam -1 . (d) First moment map of SIS-1. 
(e) Total intensity map of SIS-2. Contours are at 15.0 X (1.3, 3, 4.5, 6, 7.5, 9, 10) Jyms-^eam -1 . (f) First moment map of SIS-2. 
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(Fig. with more HI near the close companion SIS-2. All 
of the detected companions show peculiarities in their HI 
morphology. The channel maps of KUG 1359+326 reveal a 
large extension that could be tidal in origin (Figs. IS1 I10fa + 
and the galaxy's velocity field is highly distorted (Fig. \10b ) . 
Moment maps of SIS-1 (Figs. Il()b . HOt O show a kinemati- 
cally distinct feature roughly aligned with the galaxy's ma- 
jor axis, and the HI peak in SIS-2 is offset from the optical 
centre (Fig. EH 3 )- Indeed, the better we resolve a compan- 
ion to NGC 5433, the more distorted its HI morphology 
and kinematics appear. This leads us to suspect that SIS-2 
may indeed be strongly interacting with nearby NGC 5433, 
but that evidence of the encounter is concealed by our poor 
spatial and spectral resolution of this companion. 



4.3 Underlying mechanisms 

Since NGC 5433 is a starburst (see fig. 2 of Irwin et al. 1999) 
it seems plausible that the observed high-latitude extensions 
are generated by an internal mechanism. A number of sce- 
narios are consistent with the data. Symmetries in features 
about a galaxy midplane have previously been attributed to 
a single event in the underlying disc (Chaves & Irwin 2001). 
If this is the case for F2 and F3 (Fig. 0Ji), it is puzzling 
that we do not find extensions corresponding to F2 in the 
same channels as those harbouring the shell attributed to 
F3 (Fig. 03). On the other hand, the double peak in F3 it- 
self and the presence of 20 cm continuum loops in the same 
vicinity (Fig. [Tfi) suggest a common origin for all of these 
features; their locations relative to the galaxy nucleus give 
the impression of a broad-scale nuclear outflow on the west 
side of the galaxy. Alternatively, the apparent correlation 
between the continuum and HI features is also consistent 
with a galactic fountain generated by energetic events in 
the disc (e.g. Lee et al. 2001). 

None the less, one cannot dismiss the role of NGC 5433's 
active environment in the generation of the detected extra- 
planar features. For example, observations of 3 systems with 
low mass companions by van der Hulst & Sancisi (2003) all 
show evidence for gas accretion on to the parent discs. Con- 
versely, tidal interactions may pull gas out of the midplane of 
large spirals, sometimes forming arcs and plumes with sim- 
ilar morphologies to supershells (Taylor & Wang 2003). In 
the Galaxy, the origins of some high-velocity clouds (HVCs) 
have been attributed to these processes (e.g. Putman et 
al. 2003). We speculate that a similar scenario applies to 
NGC 5433: given the inferred upper limits on the masses of 
HVC analogues in external systems (e.g. Pisano et al. 2004) , 
they would not be detectable in our data. Tidal interactions 
between NGC 5433 and its smaller companions may be re- 
sponsible for the observed extraplanar emission, either via 
direct impacts, tidally induced spurs, or by fuelling the in- 
creased star formation found in the disc. Additional multi- 
wavelength observations, particularly high-resolution X-ray 
data as well as high-sensitivity H I and Ha maps, are needed 
to distinguish between these possibilities. Indeed, the envi- 
ronment can be a key component in the formation of both 
externally and internally generated high-latitude features. 



5 SUMMARY 

We have presented spatially resolved H I observations of the 
edge-on starburst galaxy NGC 5433 and its environment. 
We find that NGC 5433 has an asymmetric H I distribution 
with more gas in the southern half of the disc, and global 
properties typical of its new Sc classification. The main disc 
emission is best modelled with a gaussian ring at R = 12. 6+ 
0.3 kpc in the plane with an exponential scale-height H e ^ 
1.2 ± 0.3 kpc; the addition of a gaussian radial distribution 
with a characteristic width of D c — 3.5 ±0.3 kpc marginally 
improves the fit. 

Further examination of the data and model residuals 
reveals 3 extraplanar features, that we label Fl - F3. We 
identify coherent extensions in consecutive velocity chan- 
nels corresponding Fl and F3, and find a complete loop in 
P-V space at the location of F3. Taking this as the signa- 
ture of an expanding shell the required injection energy is 
2 x 10 54 ergs, a value comparable to that found in galactic 
and other extragalactic supershells. Now detected in sys- 
tems with a variety of star formation rates and in diverse 
environments, it appears that coherent extraplanar emission 
requiring such large input energies is a common phenomenon 
in spiral galaxies. 

NGC 5433 is in a richer environment than previously 
thought. We confirm that KUG 1359+326 is indeed a phys- 
ical companion to NGC 5433, and detect two other dwarf 
galaxies with MAPS identifiers in the field that we label SIS- 
1 and SIS-2. Including the more distant companion IC 4357 
we find that NGC 5433 is the dominant member of a group 
with at least 5 members, spanning over 750 kpc in a fila- 
mentary structure. Optical and H I distortions in NGC 5433 
itself, and in particular the pathologies in the velocity fields 
of the better-resolved companions, indicate that interactions 
are occurring in the vicinity of NGC 5433, either among the 
galaxies or with hot intergalactic gas that may be present 
in the group. 

A variety of underlying mechanisms are consistent with 
the extraplanar feature morphologies in NGC 5433. The lo- 
cations of the two peaks in F3 and their apparent correspon- 
dence with 20 cm radio continuum loops in a high-resolution 
map suggest that either a large nuclear outflow or a galac- 
tic fountain is at work. However, the active environment in 
NGC 5433 may be partly responsible for the extensions, ei- 
ther via direct impacts, tidal displacement of material from 
the midplane or by enhancing star formation in the disc. 
We argue that environmental effects may play a role in both 
internally and externally generated high-latitude extensions 
in spiral galaxies. 
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APPENDIX A: RFI EXCISION 

This appendix describes the techniques employed to remove 
strong radio frequency interference (RFI) at 1400 MHz from 
our VLA data, and our assessment of the impact of this 
RFI excision on our results. Other data reduction details 
are given in J5| 

The application of standard flux, phase and bandpass 
calibration routines to the data revealed undulations at con- 
stant declination in the central channels of all our (indepen- 
dently obtained) datasets. The culprit was the 1400 MHz 
band centre adopted for the observations: with this setup, 
correlated noise from the 7th harmonic of the VLA's 200 
MHz local oscillator at L-band contaminated the shortest 
east- west baselines in the central channels, for which the ex- 
pected fringe rate is zero (Bagri 1996). For each calibrated 
dataset, this RFI was excised by a) removing baselines with 
east-west projections close to the array centre in all chan- 
nels, and b) clipping contaminated data on amplitude in the 
visibility domain from the infected central channels only. We 
apply a) to all channels because existing AIPS software did 
not allow channel selection, and because no flux is lost from 
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the initially RFI-free channels with this approach (see be- 
low). 

The data were edited iteratively until a combination of 
parameters in methods a) and b) yielded the lowest RMS 
level a in the central channels while eliminating the RFI at 
that level. In the end, baselines with east-west projections 
within 20 metres of the array centre were flagged in each 
cube with the AIPS task UVNOU, and the optimum central 
channel clip levels for method b) were found to be 1.5 and 0.8 
Jy for the C and D configuration datasets, respectively. For 
both methods a) and b) combined, 9 per cent of the baselines 
were removed in the central channel (V = 4352 km s" 1 ) for 
the C data, and in the D data 7 per cent of the baselines were 
excised at this velocity. For all datasets, we measure a ~ 10 
per cent increase in a at V = 4352 km s -1 relative to the 
corresponding datacube mean. The same editing scheme was 
also applied in the channels adjacent to V = 4352 km s -1 
that exhibited low-level striping; in all cases less than 3 per 
cent of baselines were removed, and less than a 5% increase 
in a was measured in the resulting maps at these velocities. 

We verified that editing the RFI-free channels using 
method a) had no effect on the HI emission detected at 
those frequencies by imaging calibrated cubes in which the 
editing had and had not been applied: we find no coher- 
ent structures in a difference map of those cubes in these 
channels. To test for a systematic decrease of the HI flux 
in NGC 5433 in the central channels due to RFI excision 
with methods a) and b), we removed the same baselines 
in adjacent channels (that didn't exhibit RFI) and imaged 
them. We examined the difference between the images ob- 
tained without this additional editing and those for which 
it had been applied, in the vicinity of NGC 5433. In these 
residual maps, there is emission near NGC 5433 with flux 
at a level ~ a of the input cubes at V = 4363 km s" 1 and 
V — 4352 km s" 1 in the C-configuration data. This implies 
that the excision of RFI has resulted in a loss of flux in 
NGC 5433 at these velocities. We note that this exercise is 
strictly valid only if the HI content in the channels used in 
the test is identical to that in the central channels. At the 
resolution of our maps, the distribution of HI is similar near 
the systemic velocity of the galaxy, as is the total flux (see 
Fig- 13- It is therefore plausible that the residual map excess 
does indeed represent the extent of flux loss in NGC 5433. 

Because the potential flux loss is on the order of the map 
noise we do not correct for missing flux in any of our maps, 
nor do we expect our conclusions drawn from an examina- 
tion of HI moment maps to change. For global parameter 
computations requiring HI flux estimates of NGC 5433 in 
these central channels, we associate an additional, system- 
atic error with the measured values, at twice the difference 
map peak flux in the vicinity of NGC 5433. This system- 
atic error is represented by the dotted portion of the total 
error bars on these channels in the global profile of Fig. [5] 
We have repeated this exercise for the detected companions 
of NGC 5433 with emission in those channels in the C+D 
data. We find no evidence for any flux loss, likely due to 
their smaller extent and distance from the phase centre. 

We conclude that removing unwanted 1400 MHz inter- 
ference from our data raises a at V = 4352 kms" 1 by ~10 
per cent relative to each datacube mean, and that in the 
C configuration data the central 2 channels may bias the 
measured flux for NGC 5433 by an amount of order a (see 



Fig-HJ- All other measurements, as well as our main conclu- 
sions, are unaffected by the RFI excision applied. 

This paper has been typeset from a TjrjX/ ETfrjX file prepared 
by the author. 
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